Міністерство освіти і науки України
Національний університет “Львівська політехніка”
Звіт
про виконання лабораторної роботи № 3
з дисципліни: „Моделювання процесів
та елементів систем керування”
на тему: „ Прискорений пошук усталених режимів
електромаґнетних елементів систем керування”
екстраполяційний -алгоритм
Виконав:
студент групи КС-43
Перевірив:
Павельчак А.Г.
Львів 2008
Мета роботи – вивчити методи прискореного пошуку усталених режимів електромаґнетних елементів систем керування, а також дослідити ці режими для заданої схеми використовуючи числові методи розв’язування систем нелінійних диференціальних рівнянь.
Завдання
Завдання для лабораторної роботи дають із зазначенням номеру схеми трансформатора та варіанту чисельних значень параметрів схеми і навантаження, які необхідно взяти з лабораторної №2.
Варіант
1, Гн-1
2, Гн-1
r1, Ом
r2, Ом
r3 Ом
3
230
220
9
9
12
Варіант
C1, мФ
C2, мФ
Um, В
RН2, Ом
RН3, Ом
3
6
0,5
290
15
14
екстраполяційний -алгоритм
2.1. Домашня підготовка до роботи
1. Вивчити метод формування варіаційних рівнянь необхідних для використання моделі чутливостей до початкових умов.
2. Вивчити алгоритм побудови моделі чутливостей до початкових умов.
3. Вивчити екстраполяційний -алгоритм і метод середніх значень для пошуку початкових умов, що задовольняють умову періодичності.
4. Скласти блок-схему розрахунку стаціонарного процесу заданої схеми заданим методом.
5. Згідно блок-схеми написати програму розрахунку стаціонарного процесу заданої схеми алгоритмічною мовою С.
2.2. Робота в лабораторії
1. Ввести в комп’ютер програму, написану мовою С згідно отриманим завданням.
2. Відлагодити програму. При необхідності скоригувати блок-схему алгоритму та програму у відповідності з виявленими логічними та синтаксичними помилками.
3. Результати розрахунку стаціонарного процесу записати до текстового файлу у вигляді стовпців t, i1, i2, i3, uC1, uC2, uC3.
4. Використовуючи програму GRAFER нарисувати залежності . Можна суміщати декілька змінних однакової розмірності на одному графіку, якщо це дозволяє масштаб.
5. Остаточні версії блок-схеми, програми та отримані результати занести у звіт з лабораторної роботи.
6. Здати звіт з лабораторної роботи.
Блок-схема
Список ідентифікаторів констант, змінних, процедур і функцій, використаних у
блок-схемі алгоритму і програмі, та їх пояснення.
X_diy,Y_diy – матриці для зберігання дійсних обчислень напруг і струмів;
Xy, Yx – матриці для зберігання проміжних обчислень струмів і напруг;
K – поточний час;
OBN(t) – фунція яка обовлює тимчасові дані проміжних матриць B, Yx на дійсні;
Prurist(Dx,Dy) – фунція яка додає до даних матриць B, Yx приріст;
void Obnov_Matrix(fi) – формує матриці B, Yx виходячи з поточних струмів, напруг…
gaus – функція яка реалізує матод Гауса
fii – функція яка повертає обернене диференціальній індуктивності трансформатора, яку визначаємо за кривою намаґнечення або
e[12][12][5] – матриця яка реалізує е-алгоритм
Текст програми
Форма
Програма
#include <vcl.h>
#include <math.h>
#include <dstring.h>
#include <stdio.h>
#include <windows.h>
#include <string.h>
#pragma hdrstop
#include "Unit1.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
double x1, x2,y1,y2,m2,h,start,finish;
double a[4];
double Matrix[4][4];
double Y[4];
double X[4];
double Fai;
double Xy[5];
double Yx[5];
double B[5][5];
double Alfa[3];
double RO[3][3];
double i1,m1;
double X_diy[5],Y_diy[5];
double B_diy[5][5];
double hX[5];
double k; //
double kk; //
double C1, C2, Rn1, Rn2, Um;
double tstsrt,delt, tn, del_t_fail;
void OBN(double t);
void Prurist(double Dx, double Dy);
//---------------------------------------------------------------------------
void gaus(void)
{
double V[4][4], C[4][4],P[4];
int i,j,k;
double DY;
for(i=0;i<4;i++)
{
P[i]=Y[i];
for (j=0;j<4;j++)
V[i][j]=Matrix[i][j];
}
for (k=0;k<4...
